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This paper presents the experimental set-up of identification of harmonic 
source location in the power distribution network using time-frequency 
analysis, known as S-transform (ST) at the point of common coupling 
(PCC). S-transform offers high frequency resolution in analyzing the low 
frequency component and able to represent signal parameters in time- 
frequency representation (TFR) such as TFR impedance (Zrrr). The 
proposed method is based on IEEE Std. 1459-2010, ST, and the significant 
relationship of spectral impedances components (Zs) that been extracted 
from the Zrrr, consist of the fundamental impedance (Z7) and harmonic 
impedance (Z;). This experiment was conducted out on an IEEE 4-bus test 
feeder with a harmonic producing load in numerous different scenarios. The 
experimental was tested and verified for three consecutive months. The 
findings of this study reveal that the proposed method provides 100 percent 
correct identification of harmonic source location. 
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1. INTRODUCTION 

Due to the rise of harmonic-producing loads, harmonic distortion has become one of the primary 
power quality concerns [1], [2]. Harmonic distortion, which causes voltage and current waveforms to be 
distorted and contain various harmonic orders, is one of the most common types of disturbances [3]—[6]. The 
power system is impacted by the disturbances; therefore, monitoring is required to limit the impacts, which 
may include equipment failures due to overheating, reduced transformer life expectancy due to deterioration 
of insulation levels, nuisance tripping, and increased equipment power losses [7]—[10]. Moreover, harmonics 
can cause overheating and damage to end-user equipment, as well as have unfavorable effects on the power 
system. As a result, it is critical for a power system operator to understand the system's harmonic behavior 
[11]-[13]. As mentioned in [14]-[18], harmonic sources, on the other hand, have complicated properties such 
as nonlinearity and abrupt variations that are difficult to forecast using standard methods the foremost 
common circumstance that needs harmonic source location is to settle the disputes over who is responsible 
for harmonic distortions, whether it comes upstream or downstream of the point of common coupling (PCC) 
[19]-[21]. 
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The power direction method is the most popular method of identifying harmonic sources [22]-[24]. 
Next, the critical impedance method based on reactive power [25], [26] also offers a certain level of 
accuracy. Some basic assumptions are required for the approaches listed above, such as prior knowledge of 
source impedance [11], [27]-[29]. In contrast of active power flow direction, the reactive power methods 
provide always correct claims with regards to the dominant equivalent harmonic source. These approaches, 
however, is unable to establish the harmonic contribution of each side [30]—[34]. Other methods, such as 
fluctuation and regression methods, require that the major harmonic source be on the customer side and that 
the background harmonic voltage required to be stable [35]-[39]. Furthermore, approaches based on the 
detection of total harmonic distortion are insensitive to changes in the phase angle of the harmonic source, 
making it impossible to precisely establish the source of harmonic distortion [40]-[42]. According to [43], 
current techniques have been used to identify the harmonic contributions of the customer and the utility in 
order to detect the harmonic source. Based on the reference impedance as in [34], [44]-[47], a harmonic 
vector approach has been suggested to determine the utility and customer's harmonic contributions at the 
PCC. This method allowed for the calculation of harmonic contributions without determining customer 
impedance, and also improved the findings in resonance situations. The independent component analysis 
(ICA) methods were utilized in recent research [47]-[50], which need the impedance on the customer's side 
to be higher than the one on the utility side. However, when the network contains filters or capacitors on the 
customer side, this is impractical [7], [51]. 

Short-time fourier transform (STFT) and Stockwell transform or S-transform (ST), are the most 
common time-frequency domain transforms employed in harmonic signal detection approaches [52]-[58]. 
As explained in [59]—[62], STFT on the one hand, has some disadvantages such as, this transform is window- 
dependent and has a fixed resolution based on the window size. Furthermore, because the STFT is an Fourier 
transform (FT) based technique, it may have issues with the picket-fence effect [63], [64]. ST, on the other 
hand, because it is a multiresolution spectrum analysis technique, does not have these issues [65]—[68]. As a 
result, ST appears to be a potential transform for power system protection [27]. Because it incorporates 
information in both the temporal and frequency domains, the ST has shown to be effective in harmonic signal 
identification approaches [69]. Thus, this paper proposes an experimental setup based on IEEE Std.1459- 
2010 and ST due to identify the harmonic source location. 


2. METHOD 
2.1. Proposed method 

The identification of harmonic source location is divided into five steps, as indicated in Figure 1. 
The signals are first measured for both voltage and current at the PCC of the network system. Second, four 
specific instances were explored for recognising harmonic sources on IEEE 4-bus test feeders [44]. The time- 
frequency representation (TFR) analysis was done on the PCC's voltage and current measurements in the 
third step (Vpcc and Ipcc). This analysis yielded the impedance TFR (Zrrr), which was then used to calculate 
the impedance spectral (Zs) components by calculating the values of the Zrrr components. Finally, the 
significant association between the fundamental impedance (Z;) and harmonic impedance (Z;,) components of 
impedance spectrum (Zs) components was observed and employed for harmonic source detection. In this 
experiment, a harmonic generating load was chosen with an amplitude modulation index (ma) of 1.0, a 
frequency modulation index (ms) of 90, and an input frequency (fi) of 50 Hz [70]-[73]. 


Voltage and Current Measurement at 
the PCC 


N-N N-H H-H H-N 
—] 
$ 4 
Voce Ipce 
Voltage TFR Current TFR 
t ] t 
Impedance TFR 


¥ 


Harmonic Source Identification 


Figure 1. The implementation of the harmonic source identification method 
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The IEFE 4-bus test feeder is chosen and illustrated in Figures 2 and 3 in order to detect the 
harmonic source location in the power distribution network in consideration of upstream and downstream of 
the PCC. Where N is a non-harmonic source which is the resistor load and H is the harmonic producing load. 
In order to test and evaluate the proposed method, an experimental setup was built up in an advanced digital 
signal processing (ADSP) research facility, as illustrated in Figure 4. 


Figure 4. Experimental setup of the proposed system 


2.2. S-transform 

The S-transform (ST) is hybrid of wavelet transform (WT) and the STFT, which inherits the 
advantages of both in signal processing [74], [75]. In the transformation process, ST uses a moving and 
scalable localising Gaussian window in particular. ST can be defined as shown in: 


-(t-t)2 f2 
2 


STEA =S x() e ee (1) 
1 1 =e? 
o(f) = 3 g(t) = Speer? (2) 


where x(t) is the signal, t is the time, f is the frequency, g(t) is the scalable Gaussian window, and o(t) is a 
parameter that controls the position of the Gaussian window. 


2.2.1. Signal parameters 

The TFR is used to determine the signal parameters of power quality. Furthermore, the 
instantaneous value is used in the analysis to obtain real-time parameters [76]. 
— Instantaneous root-mean square voltage 

The root-mean square (RMS) voltage of signal (Vins(t)) can be obtained from the sampled 
waveform, and written as [77], [78], 
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Vms (t) = IE P.(t, f) af (3) 


— Instantaneous root-mean square fundamental voltage 
The instantaneous RMS fundamental voltage (Virms(t)) can be computed as [79], 


Varms(t) = [2 J7™ Pet, f) df (4) 


fni = fo + 25 HZ; fio = fo — 25 Hz 


where Px is the power spectrum obtained from the TFR of signal and fo is the fundamental frequency 
corresponding to the power system frequency. 


2.2.2. Impedance time-frequency representation analysis 

The impedance TFR (Zrrr) offered useful information about the frequency response of the system, 
as well as harmonic points and possible issues caused by harmonic distortions. The desired current harmonic 
data and the difference in voltage harmonic data at the location of interest have to be measured in order to 
determine the impedance TFR. The Zrrr at each harmonic frequency was calculated using this data, and the 
results were shown [80]. The Zrrr equation can be expressed as (5), 


Svt f) 
Z 5 
TFR Sit f) ( ) 


where Sy(t f) signifies the TFR of voltage and S,(t,f) signifies the TFR of current. 


2.2.3. Spectral impedance 

The spectral impedance (Zs) comprises the fundamental impedance (Z;) and harmonic impedance 
(Zn) that is obtained from Zrrr [81]. The fundamental impedance (Z;) was an impedance at 50 Hz, which was 
the frequency of the power source. In the meantime, harmonic impedance (Z) was a harmonic impedance 
with an order of harmonics. 


3. RESULTS AND DISCUSSION 

The implementation of the proposed technique initially done by measuring the voltage and current 
signals at PCC with consideration of 4 specific cases as discussed in 2.1. The linear time-frequency 
distribution method namely S-transform is applied in the analysis. The location of harmonic sources can be 
distinguished by analyzing the significant relationship between Z; and Zn, accordingly. 


3.1. Case 1: No harmonic source 

Only the linear loads were placed upstream and downstream of the PCC in case 1. The voltage 
signal in the time domain, as well as its voltage TFR, are shown in Figures 5(a) and 5(b). The maximum 
voltage was 342.5 V, while the maximum current was 66.5 A. Meanwhile, Figures 5(c) and 5(d) illustrate the 
TFR of voltage and current signals derived from S-transform analysis. The higher the magnitude, the redder 
the colour bar, the lower the magnitude, and the bluer the colour bar. There were no other components in the 
signals and the largest magnitude was only seen at 50 Hz. The results showed that there were no harmonic 
components in the signal. The Z; existed at 50 Hz with a resistance of 4.8 ohm and no harmonic components 
in the signal, as shown in Figure 5(d). Thus, in case 1, the significant relationship between Z; and Z in the 
power system network at no harmonic producing load can be expressed as, 


Zi + 0 ohm (6) 
Zn =0 ohm (7) 


where for harmonic component, h is any positive integer, whereas for interharmonic, h is any positive non- 
integer. 
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Figure 5. Case 1: (a) voltage signal in time domain, (b) current signal in time-domain, (c) TFR impedance 
using S-transform, and (d) spectral impedance 


3.2. Case 2: Harmonic source located at point of common coupling’s downstream 

The linear load is positioned upstream of the PCC in case 2, while the harmonic load is located 
downstream. The TFR of voltage and current signals derived from the S-transform analysis is shown in 
Figure 6(a) and (b). It can be seen that the harmonic and interharmonic components exist between 200 and 
1000 Hz, whereas the fundamental components of voltage and current have the maximum magnitudes at 50 
Hz. The Zrrr is calculated using (5) in Figure 6(c), and the figure demonstrates that impedance components 
occur at frequencies of 50 Hz, 275 Hz, 375 Hz, 600 Hz, 700 Hz, and 900 Hz, respectively. The Zs is then 
derived by calculating the parameters of the Zrrr, as shown in Figure 6(d). 

Table 1 summarises the Zs characteristics shown in Figure 6(d). The Z; value is always higher than 
any Zn components, as can be seen. The relationship between the Zs components can be used to identify the 
location of harmonic sources, according to the findings. As a result, in instance 2, the significant relationship 
between Z; and Z, at the condition of the harmonic source downstream of the PCC can be stated as, 


Zi + 0 ohm (8) 
Ln < Zi (9) 
where for harmonic component, h is any positive integer whereas, for interharmonic, h is any positive non- 


integer. 


Table 1. The spectral impedance components for case 2 
Spectral impedance | Ohm 


Zı 3.4 
Z275 25 
Z375 2.0 
Z600 2.1 
Zr00 1.5 
Zo00 2.0 
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Figure 6. Case 2: (a) voltage signal in TFR using S-transform, (b) current signal in TFR using S-transform, 
(c) TFR impedance using S-transform, and (d) spectral impedance 
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3.3. Case 3: Harmonic sources located at point of common coupling’s upstream and downstream 

The TFR of voltage and current signals derived from the S-transform analysis for case 3 is shown in 
Figures 7(a) and 7(b). It can be seen that the harmonic and interharmonic components occur in the frequency range 
of 200 Hz to 1000 Hz, whereas the fundamental components of voltage and current, respectively, have the 
maximum magnitudes at 50 Hz. The voltage and current waveforms can be seen to be distorted due to the harmonic 
load located upstream and downstream of the PCC. The Zrre is calculated using (5) in Figure 7(c), and the figure 
demonstrates that impedance components occur at frequencies of 50 Hz, 275 Hz, 375 Hz, 600 Hz, 700 Hz, and 900 
Hz, respectively. The Zs is then calculated by estimating the parameters of the Zrrr, as shown in Figure 7(d). 
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Figure 7. Case 3: (a) voltage signal in TFR using S-transform, (b) current signal in TFR using S-transform, 
(c) TFR impedance using S-transform, and (d) spectral impedance 
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Table 2 summarises the ZS characteristics shown in Figure 7(d). The Z1 value is always lower than 
any Zh components, as can be shown. The relationship between the Zs components can be utilised to pinpoint 
the location of harmonic sources, according to the findings. As a result, in case 3, the significant relationship 
between Z; and Z, at the condition of the harmonic source positioned upstream and downstream of the PCC 
may be expressed as shown in: 


Zi + 0 ohm (10) 
Ln > Zi (11) 
where for harmonic component, h is any positive integer whereas, for interharmonic, h is any positive non- 


integer. 


Table 2. The spectral impedance components for case 3 
Spectral impedance Ohm 


Zı 0.7 
Z275 1.6 
Z375 1.8 
Z600 1.9 
Z700 1.4 
Zo00 2.1 


3.4. Case 4: Harmonic source located at point of common coupling’s upstream 

The harmonic load is positioned upstream of the PCC in case 4. Figures 8(a) and (b) illustrate the 
voltage and current signals acquired from S-transform analysis in the TFR. Between 200 Hz and 1000 Hz, 
the lowest-magnitude harmonic and interharmonic components are present, with the maximum component 
magnitude at 50 Hz. The Zrrr is calculated using equation 5 in Figure 8(c), and the figure demonstrates that 
impedance components occur at frequencies of 50 Hz, 275 Hz, 375 Hz, 600 Hz, 700 Hz, and 900 Hz, 
respectively. The Zs is then calculated by estimating the parameters from the Z7rr, as shown in Figure 8(d). 
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Figure 8. Case 4: (a) voltage signal in TFR using S-transform, (b) current signal in TFR using S-transform, 
(c) TFR impedance using S-transform, and (d) spectral impedance (Zs) 
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Table 3 summarizes the Zs characteristics shown in Figure 8(d). The Z; value is the same for all Z, 
components, as can be observed. At the condition of the harmonic source positioned upstream of the PCC, 
the significant relationship between Z; and Z; can be stated as, 


Zi +0 ohm (12) 
Zr-Zı (13) 


where for harmonic component, h is any positive integer whereas, for interharmonic, h is any positive non-integer. 


Table 3. The spectral impedance components for case 4 
Spectral impedance Ohm 


Zi 3.3 
Z275 3.3 
Z375 3.3 
Z600 3.3 
Z700 3.3 
Zo00 3.3 


Furthermore, the proposed method was tried and verified on an experimental setup in October, 
November, and December 2021, with the harmonic producing load in the linear area (amplitude modulation 
index is 0 < ma < 1 and inverter switching frequency range is between 2 kHz and 15 kHz) [82]. Surprisingly, 
as demonstrated in Figure 9, the proposed method offers 100 percent accurate harmonic source location 
detection. According to Table 4, the proposed method is 100 percent correct in each scenario, and the 
significant relationship of Zs for harmonic source location identification is summarized as shown in: 


Percentage 


Figure 9. The correctness of the proposed method 


Table 4. Result of proposed method 
Z, Zn Remark 
Z;#0ohm Z,=0 Case 1: No harmonic source 
Z,#0ohm Z, < Z, Case 2: Harmonic source located at PCC downstream 
Z;#0ohm Z,>Z; Case 3: Harmonic source located at upstream and downstream of PCC 
Z,;#0ohm Z,=Z, Case 4: Harmonic source located at PCC downstream 


4. CONCLUSION 

Time-frequency distribution analysis namely S-transform has shown tremendous result in this 
analysis. The major contribution of this study is the discovery of a significant relationship between Zs 
components that acquired from S-transform analysis in locating harmonic source location. As a result of the 
proposed method's results, the harmonic source site can be identified using the significant relationship of 
spectral impedances in a fast, cost-effective, and accurate manner. 
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